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Abstract 

We discuss corrections to the Leading-LogQ 2 relationships between the gluon density g(x,Q 2 ) 
and F L (x,Q 2 ), dF T {x,Q 2 )/dlogQ 2 in the HERA range of large ~. We find that the above 
quantities probe the gluon density g(x, Qtl) a ^ Qtl = Ct,lQ 2 , with the Q 2 -rescaling factors 
Ct ~ 2.2 and Cl ~ 0.42. The possibility of treating charm as an active flavour is critically 
re-examined. 
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1 Introduction. 



At large - , the gluon desnity g(x, Q 2 ) is much higher than the density of charged partons 
q(x, Q 2 ), q(x, Q 2 ) (hereafter x is the Bjorken variable and Q 2 is the virtuality of the photon) and 
the photoabsorption will be dominated by interaction with the nucleon of the lightcone qq Fock 
states of the photon via the exchange by gluons (Fig. 1) or, alternatively, by the photon-gluon 
fusion 7*g — ► qq. Consequently, the longitudinal structure function F L (x,Q 2 ) and the slope of 
the transverse structure function Ft(x,Q 2 ) = 2xFi(x,Q 2 ) become direct probes of the gluon 
structure function G(x, Q 2 ) = xg(x, Q 2 ): 

F L (x,Q 2 ) = ^pY.e}G(Z L x,Q 2 ), (1) 

In Eqs. dDH) OisiQ 2 ) is the running strong coupling, e/ is the quark charge in units of the 
electron charge and the summation goes over the active flavours /. The suggestion of Fl as a 
probe (partonometer) of the gluon density is due to Dokshitzer [1], Eq. ([J) was elaborated in 
[2]. Eq. (p|) readily follows from formulas (6.37) and (6.34) of the Roberts' textbook [3] (see also 
[4]). Both equations were derived in the Leading-LogQ 2 approximation (LLQA), the x-rescaling 
factors £t,l ~ 2 [2-4]. The emergence of the gluon-dominated scaling violations at large - was 
clearly demonstrated in the recent QCD-evolution analysis of the NMC structure functions [5]. 
The obvious advantage of Eqs. (0,0) is that one does not need solving the coupled QCD-evolution 
equations for the gluon and (anti) quark densities. 

As a byproduct of our analysis [6] of determination of the BFKL pomerons intercept from 
Fl(x, Q 2 ),d Ft (x,Q 2 )/d log Q 2 we have noticed the potentially important corrections to the 
LLQA. One obvious issue is which Q 2 is sufficiently large for charm to be treated as an active 
flavor, because for the (u, d, s) active flavours e / = § compared to J2 e / — y if charm also were 
an active flavour. Our analysis [6] suggests that it is premature to speak of Nf = 4 active flavours 
unless Q 2 > (100 - 200)GeV 2 . In [6] we also noticed that F L (x,Q 2 ) and dF T (x, Q 2 )/d\ogQ 2 
probe the gluon structure function G(x, Qtl) a ^ different values of Q\ L , which are both different 
from Q 2 . Because of importance of determination of the gluon structure function, which is a 
fundamental quantity in the QCD parton model, in this paper we present an update of Eqs. ([1],|2]) 



which are valid also in the BFKL (Balitzkii-Fadin-Kuraev-Lipatov [7]) regime, i.e., beyond the 
LLQA. 

At large ^, deep inelastic scattering can be viewed as an interaction with the target nucleon 
of the lightcone qq Fock states of the photon via the exchange by gluons (Fig.l). The principal 
quantities are the probability densities \^t,l{z, r)\ 2 for the qq Fock states with the transverse size 
r and the fraction z of photon's lightcome momentum carried by the (anti)quark, and a(x,r) - 
the total cross section of interaction of the qq colour dipole of transverse size r with the nucleon 
target [8,9]. This dipole cross section satisfies the generalized BFKL equation derived in [9,10] 
and is related to the differential density of gluons by the equation [9,11,12] 



7i , . , r d 2 k 4[1 - exp(ikr)] dG(x a , k 2 ) , . 

The wave functions of the qq Fock states of the (T) transverse and (L) longitudinal photon were 
derived in [8] and read 



\* T (z,r)\ 2 = EeJI^'W)! 2 = *^£ e J{[*» + (1 " zWMer? + m 2 f K (er) 2 } , (4) 



\* L (z,r)\ 2 = £4l4 /7 W)| 2 = £^J>} Q 2 z 2 (l - z) 2 K (er) 2 , (5) 



(2tt) 2 \ 

where K v (x) are the modified Bessel functions, e 2 = z(l — z)Q 2 +m 2 and rrif is the quark mass. 
The resulting photoabsorption cross sections are equal to ([8], see also [13]) 

ar^N^Q 2 ) = f Q dz J d 2 r\* T (z,r)\ 2 a(x,r) = ^£e 2 ^ / ^ 



/ 



r d 2 K I [z 2 + (1 - z) 2 ]k 2 + m 2 [z 2 + (1 - z) 2 ]k(k + k) + m? f \ dG{x g , k 2 ) 2 
1 K 2 + e 2 \ k 2 + e 2 {k + nf + e 2 J d\ogk 2 a (6) 

a L ( 7 *iV,x,Q 2 ) = f Q dz J d 2 r\^ L (z,r)\ 2 a(x,r) = ^ £ e) £ dzAQ 2 z 2 {\ - zf J £± 

r d 2 n i 1 1 )dG(x g ,k 2 ) 2 

7 ^T^lF + ^~jfc(jfc + «)+e2j Slog A; 2 U 

Here the running coupling 0:5 (g 2 ) enters the integrand at the largest relevant virtuality, 

q 2 = max{£ 2 + k 2 , k 2 } , (8) 
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and the density of gluons enters at 



x g = x(l + ^) (9) 
where M t is the transverse mass of the produced qq pair in the photon-gluon fusion 7*g — > qq: 



M? 



„ m 2 f + k 2 m 2 f + (k + k) 2 , 
M 2 = s + -L — L '- . (10) 



The flavour and Q 2 dependence of structure functions is concentrated in wave functions 
whereas the dipole cross section a(x,r) (the differential gluon density dG(x,k 2 )/d\ogk 2 in the 
momentum representation) is universal for all flavours. The important virtue of the (r, z) rep- 
resentation in (fU?]) is the factorization of integrands as \*&t,l(z, r)\ 2 a(x, r), which corresponds 
to the exact diagonalization of the diffraction scattering matrix in the (r, ^-representation. 
Furthermore, the above dipole-cross section representation (||[7p and wave functions (|,|]) are 
valid also in the BFKL regime, i.e., beyond the LLQA, and allow an easy incorporation of 
the unitarity (absorption) corrections at large - [9,11], with allowance for which Eq. (^) must 
rather be regarded as an operational definition of the gluon density beyond the LLQA. The x 
(energy) dependence of the dipole cross section a(x,r) comes from the higher qqg\--.g n Fock 
states of the photon, i.e., from the QCD evolution effects, described at large - by the gen- 
eralized BFKL equation [9,10]. The structure functions are given by the familar equation 
Ft,l(x,Q 2 ) = (Q 2 /4:7ia em )(jT,L(x,Q 2 )- We advocate using Fl and Ft = 2xFi, which have 
simpler theoretical interpretation compared to F 2 = Fx + F^ which mixes interactions of the 
transverse and longitudinal photons. For the sake of completeness, we notice that, in the Born 
approximation the differential gluon density is related to the two-body formfactor of the nucleon 
(N\ exp(iA; 1 r 1 + ik 2 r2)\N) by the equation [12] 

^ (X 'f 2 2) = -a s (k 2 )(l - (N\eMik(n - f 2 )\N) , (11) 

O log k z 7T 

and the limiting form of Eqs. (||[7|) derived in [8] is obtained. 

The leading contribution comes from values of M 2 ~ Q 2 , so that x g ~ 2x. Strictly speaking, 
the x-rescaling factors can not be determined within the Leading-Log^ approximation. From the 
practical point of view, when analysing Fl and 8Ft / d log Q 2 , it is sufficient to use x g = £,t,l% 
with C,x,l ~ 2 as determined in [2-4]. In this communication we concentrate on corrections to 
the LLQA. 



2 Active flavours and the onset of LLQA 

The ratio a(x ) r)/r 2 is a smooth function of r. Similarly, dG(x, k 2 )/d\ogk 2 is a smooth function 
of k 2 . Consequently, it is convenient to use the representations 



1 r dr 2 <j(x, r 



FtM 2 ) = -J ^-^^W\Q 2 y) ^ (12) 



ir 3 J r 2 r 2 ^ f 



= ^E4/feWV,F)^!), (i3) 



^^ = i/^^E^«V) (i4) 

where the weight functions $^ ^ and Wjfp are defined by 

$V f ~\Q 2 ,r 2 ) = (7r 2 /4a em ) C dzQ 2 r A \^ J \z,r)\ 2 , (16) 

Jo 

wf /_) (Q 2 ,r 2 ) = (vr 2 /4a em ) C dz Q 2 r*\¥/ f ~\z,r)\ 2 , (17) 

Jo 

9 log Q 2 



^ (QVH ?£pM (18) 



and the kernels 0t,l are given by 

r 2 n. v (Q 2 ) :>r'H 



e T ,,(Q 2 ,* 2 ) = gV2)> 



Here Jo(x) is the Bessel function and the running coupling as {q 2 ) enters at the largest relevant 
virtuality: q 2 = max{A; 2 , C 2 /r 2 }, where C ~ 1.5 ensures the numerically similar results of 
calculations in the (r, z) and the momentum representations [8] . Since the BFKL equation is 
only known to the leading order, and the differences between the leading-order and next-to- 
leading order at large Q 2 are marginal [3], below we use the one-loop strong coupling as{k 2 ) = 
4-7T / (3 \og(k 2 j A 2 ) with A = 0.3GeV. Here (3 = 11 — \Nf, and in the numerical estimates we 
impose the infrared freezing as{k 2 ) < = 0.8. 
The conventional LLQA corresponds to 

4 //) (Q 2 ,r 2 )=^(Q 2 -i), (20) 



Q« f L \Q 2 ,k 2 )=0(Q 2 -k 2 ) (21) 

and to the factoring out as{q 2 ) ~ &s{Q 2 ) from the integral (|19|). A derivation of the LLQA 
formulae (JUID is based upon precisely these approximations. Our analysis of approximation (|20|) 
in [6] revealed a very slow onset of LLQA for the charmed quarks. Below we analyse in more 
detail an accuracy of Eq. (|2lD and of the LLQA relations 

Our results for Qtl(Q 2 > are shown in Fig. 2. We find a very strong departure from the 
LLQA Eq. fl2"T|): i) the kernels ©t,l have a very broad diffuse edge, ii) the position of the diffuse 
edge is shifted compared to the naive expectation k 2 = Q 2 , iii) the kernels Qt,l flatten at large 
Q 2 , but the height of the plateau is different from unity, iv) the onset of LLQA for charm is very 
slow. Discussion of these effects is particulary simple in terms of the representation (|T9|). 

In order to facilitate the further discussion, we remind the salient features of weight functions 
W^ f P{Q 2 ,r 2 ) [6]. Firstly, at asymptotically large Q 2 , 

J^ r W«P(Q 2 ,r 2 ) = l. (22) 

Secondly, W^PiQ 2 , r 2 ) are peaked at r 2 = B TtL /Q 2 , where 

B T ^2.3, (23) 
B L ^11. (24) 

Therefore, F L (x,Q 2 ) and dF T {x,Q 2 )/d\ogQ 2 probe a(x, r) at r 2 = B L)T /Q 2 - Notice a sub- 
stantial departure from the naive LLQA expectation of Bt,l ~ 1- The width of these peaks is 
quite broad, A\og(Q 2 r 2 ) ~ 3 [6]. Furthermore, the function f(x) = 4[1 — J (x)]/x 2 , shown in 
Fig. 3, is similar to the step-function, but has a very broad diffuse edge. The position of the 
edge corresponds to an effective step-function approximation f(x) = 6(A a — x) with A a w 10, 
and leads to a dramatic deviation from the naive LLQA estimate Q 2 ~ 1/r 2 in the small-r limit 
ofEq. (|): 



vr 2 , , „, „ A 
3 ^ — v^-v~'- r 2 

Now we discuss the results for Qt,l{Q 2 , k 2 ) in the light of these observations. 



The diffuse edge of kernels ©t,l originates from the two factors: a relatively large width of 



the weight functions Wt,l(Q 2 ,^ 2 ), and a broad diffuse edge of f(kr) = 4[1 — J (kr)]/(k 



r) 2 . 



The kernels ©t,l flatten at k 2 <C Q 2 , but the height of the plateau H^l L (Q 2 ) is significantly 
different from unity. At small Q 2 the height of the plateau is smaller than unity, which signals 
the sub-LLQA for the considered flavour. This effect is particularly important for heavy flavours 
(charm, bottom). For instance, for the charmed quarks Hrf^ only very slowly rises with Q 2 
reaching only ~ .5 at Q 2 = 30GeV 2 , so that charm is an only marginally active flavor and 
for the charm conrtribution to the structure function the LLQA is badly broken unless Q 2 > 
(100 - 200)GeV 2 . For the 6-quarks the height of the plateau is « 0.3 at Q 2 = 120GeV 2 and 
~ 0.7 at Q 2 = 480GeV 2 . The onset of the LLQA for heavy flavours is particularly slow in the 
case of the longitudinal structure function, because merely by gauge invariance the longitudinal 
cross section is suppressed, (Tl/&t oc Q 2 /Am 2 at Q 2 <4mj (for instance, see [14]). 

On the other hand, the excess over unity comes from as{q 2 ) / as{Q 2 ) > 1 in the integrand of 
Eq. flTjp. Of course, at the asymptotically large Q 2 , the height of the plateau H^[/{Q 2 ) tends 
to unity, because here the LLQA becomes accurate and one would have replaced as(<f) i n the 
integrand by as{\) ~ OisiQ 2 )- Closer inspection of Eq. (|T9| ) shows that, at the moderately 
large Q 2 , besides the width of Wt,l(Q 2 \ r2 ) and the diffuse edge of f(x), a large contributor to 
as(q 2 ) > OisiQ 2 ) is the large value of A a . In the region of Q 2 of the practical interest for the 
HERA experiments, the excess of H^[\Q 2 ) over unity is particularly large for the light quarks. 

At sufficiently high Q 2 > 4m 2 , the position of the diffuse edge of &t,l(Q 2 , k 2 ) in the natural 
variable k 2 /Q 2 is approximately flavour- independent. Notice, that the diffuse edge is definitely 
shifted towards positive \og(k 2 /Q 2 ) = log(CV) ~ 1 for the slope of the transverse structure 
function, and for the longitudinal structure function the position of the edge is shifted towards 
negative \og(k 2 /Q 2 ) = log(Cy ~ —1. We find strong departure from the LLQA assumption 
Ct,l = 1- 

3 Measuring the gluon distribution 

The differential gluon structure function dG(x,Q 2 )/d\ogQ 2 is a slow and smooth function of 
logQ 2 . For this reason, we can quantify the shift of the diffuse edge and the height of the 
plateau of global kernels ®r,L(Q 2 ,k 2 ) for Nf = 5 flavours (u,d, s,c,b) in terms of the effective 



7 



step-function parameterization 



®tAQ\ fc2 ) = 77 E epT,i(Q 2 , k 2 ) = ^-H T , L {Q 2 )d{C T , L {Q 2 )Q 2 - k 2 ) . (26) 

I 1 /=1 L 



The so-determined height Ht,l{Q 2 ) of the effective step-function and the Q 2 -rescaling factor 
Ct,l{Q 2 ) are shown in Fig. 4. They enter the modified relationships (0]|^) as follows: 

F L (x, Q 2 ) = ^l^H L (Q 2 )G^ L x, C L {Q 2 )Q 2 ) (27) 
dF T (x,Q 2 ) _ a s (Q 2 ) 11 (n ^ r( t ^ r (n 2. n 2, (9R , 

d\o g m ~ ~^t^ Ht{q )g ^ tX ' Ct{q )q ] (28) 

The viable approximations for the kinematical range of the HERA experiments are Ct{Q 2 ) = 2.2 
and Cl(Q 2 ) = 0.42. Because of variations of the factor a.s{q 2 ) / ots{Q 2 ) in the integrand of (|T9|). 
these results for Ct,l slightly differ from, but are close to, the estimate Ct,l{Q 2 ) ~ Bt,li 'A a . 

In the case of the slope of the transverse structure function, there is a significant numerical 
cancellation of effects of the suppressed, sub-LLQA contribution of heavy flavours, and of the 
enhanced height of the plateau for light flavours. Because of this numerical conspiracy, despite 
the charm not being an active flavour at all, we find Ht(Q 2 > lOGeV 2 ) = 1 to better than 10 per 
cent accuracy in the kinematical range of HERA. For this reason, the recent estimate [15] of the 
gluon structure function from the HI data on the scaling violation, using the the Nf = 4 LLQA 
formula (^|), must be regarded as numerically reliable. For the longitudinal structure function 
this cancellation is less complete. The residual departure of Hl(Q 2 ) from unity is a monotonic 
function of Q 2 and remains quite substantial in the kinematical range of HERA . 

The accuracy of relations (|27|j28| ) can be checked computing first the structure functions for 
certain parametrization of the gluon density, and then comparing the input gluon density with 



the ouptut from ( p7| , |2~8"D . We have performed such a test for the gluon densities [11,12] and [16], 
both of which give a good description of the HERA data on the proton structure function [17]. 
For the derivative of the transverse structure function, the input/output comparison suggests 
that the accuracy of the above procedure is as good as 5 per cent at Q 2 > lOGeV 2 in the whole 
HERA range of x. For the longitudinal structure function, the input /output agreement is better 
than 10 per cent at x > 10~ 4 and Q 2 > lOGeV 2 , but gets worse at smaller x and not so large 
Q 2 . At x = 10 5 the 10 per cent agreement only holds at Q 2 > 30GeV 2 . These estimates of the 



S 



accuracy of relations ( p7| , p8| ) can easily be improved with the advent of high accuracy data on 
the gluon structure functions from the HERA experiments. Small corrections from the z — r 2 
correlations in the wavef unctions and/or the x g — (k, k) correlations implied by Eqs. ([9"|,|IUD, 
can also be easily included should the accuracy of the data require that. 



4 Conclusions 

We re-examined determination of the gluon structure function G(x, Q 2 ) from the longitudinal 
structure function Fl(x, Q 2 ) and from the slope 8Ft(x, Q 2 )/dlogQ 2 of the transverse structure 
function, and derived new relationships Eq. (^) and Eq. (|28|) . In the range of Q 2 of the interest 
for HERA experiments, charm is only marginally active flavour. None the less, because of 
the numerical conspiracy of corrections to LLQA for the light and heavy flavour contributions, 
the overall normalization factor Ht{Q 2 ) ~ 1 as if all the Nf = 5 flavours (u, d, s, c, b) were 
active. Here the major effect is the Q 2 -rescaling Ct ~ 2.2. We find that the onset of LLQA 
for the longitudinal structure function is much slower than for the transverse structure function. 
Numerical experiments suggest that our relationships ( p7| , p8D have < 10% accuracy. 
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Figure captions 

Fig.l - Leading QCD subprocesses at large K 

Fig. 2 - The kernels Qt,l(Q 2 , k 2 ) for the light quarks (u,d), the charmed quark and the global 
kernel for Nf = 5 flavours (u,d,s,c,b). The dot-dashed, dashed, solid, double-dot-dashed 
and dotted curves are for Q 2 = 0.75, 2.5, 8.5, 30, 480GeV 2 , respectively. 

Fig. 3 - The function f(x) — 4[1 — J (x)}/x 2 . 

Fig. 4 - The height of the plateau H TtL (Q 2 ) and the C} 2 -rescaling factor C TjL (Q 2 ) as a function 
of Q 2 . 
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